Magnetic oscillations and frequency mixing in a two-band conductor 
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Abstract 

Exact analytical results of the de Haas-van Alphen (dHvA) effect in an idealized two-band Fermi liquid with 
parabolic dispersion are presented. We consider a Fermi surface consisting in two electron bands with different 
band edges and band masses. Magnetic breakthrough (MB) between the bands is negligible. Analytical ex- 
pressions of the dHvA Fourier amplitudes are derived in the case where the total number of electron is fixed 
(Canonical Ensemble, CE). As already reported in the literature, the oscillations of the chemical potential 
yield frequency mixing and Lifshitz-Kosevich (LK) theory, which is valid in the Grand Canonical Ensemble 
(GCE), does not apply at very low temperature. We show that the corresponding Fourier amplitudes de- 
pend on the commensurability between the two effective masses and also the two fundamental frequencies. 
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Experimental data on dHvA quantum oscilla- 
tions in organic conductors show evidence of a devi- 
ation with LK [1] and Falikov-Stachowiak (FS) [2] 
theories for low dimensional systems. We observe 
in the Fourier spectrum of the magnetization in the 
transfer salt k - (BEDT - TTF) 2 Cu(NCS) 2 [3] 
frequencies that are not related to a classical or- 
bit on the Fermi surface. For these "forbidden" 
frequencies, the electron orbits locally in the oppo- 
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site direction imposed by the magnetic field. This 
material is made of 2 disconnected Fermi surfaces, 
one pocket hole and an open one dimensional 
surface. From LK and FS theories, we should ob- 
tain a small frequency f a corresponding to the 
pocket, and, by MB effect between the pocket and 
the open surface, a larger frequency fp > f a at 
higher field (>20T at IK), plus their harmonics 
and classical combinations that are permitted like 
fa + fp, 2fp - fa etc... due to MB and Bragg 
reflections throughout the Brillouin zone. The 
"forbidden" frequencies seen experimentally cor- 
respond to ff3 — fa, 2ff3 — 2 f a - In Shubnikov-de 
Haas (SdH) transport experiments, they may be 
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explained by quantum interference or Stark effect, 
like in {BEDO - TTF) 5 [CsHg{SCN) 4 } 2 [7]. In 
dHvA experiments, the explanation seems rather 
different because Stark effect does not apply. 
There is an argument [4,5,6] showing that fluctu- 
ations of the chemical potential are important in 
these low dimensional systems, especially at very 
low temperature (<1K) in the high field range. 
This implies that LK and FS frequencies and 
Fourier amplitudes should be different in the CE, 
where the electron density n e is fixed, and that 
the classical way of defining a frequency as area 
of a classical orbit is incorrect in this ensemble. 
However, since the LK and FS are correct in the 
GCE, thermodynamical relations connect the 2 
Ensembles. Determinating the chemical potential 
is not simple in the finite temperature case and a 
complete calculation is necessary. However we can 
obtain analytical results at T=0 in some cases [6] . 

In order to reproduce the main features of 
transfer salt Fermi surface, we condidcr a two- 
band model with gaps A Q= o.i and effective masses 
m a=01 in an external magnetic field B (we ne- 
glect magnetic interaction). This model has been 
considered by other authors mostly in a numerical 
way [4,8]. However there is no MB unlike BEDT 
systems. The Landau levels (LL) of each band 
a = 0, 1 have energies: 



e a (n) 
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We will define energies in unit of 2irh 2 /mj, and 
b = eB/h is the dimcnsionless field per unit cell 
area. In absence of field, the electron densities 
n a of each band satisfy the conservation equation 
n « — n e- If S a is the area of the ce-band Fermi 
surface, then n a = S a /Air 2 . In term of the Fermi 
energy Ef, we have S a = 2irm* a /Ti 2 (Ef — A a ). In 
the following, we will use the dimensionlcss vari- 
able x — n e /b, and the ratios c a — TOp/m* . In the 
new units, the LL are e a (n) = A a + c a b(n + 1/2). 
It is also convenient to use the dimensionless 
fundamental frequencies f a — S a /A-K 2 n e , with 
fo + fi = 1- Defining f b = (A - Ai)/cm e , we 
obtain 



fo = 
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(1 - f b ), and /i = 
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From the LK theory, the Grand Potential fl is 
the sum of each band contribution Q : 
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Ra,p(T) — p\ a j sinh(pA a ) is the temperature 
reduction factor, and X a — 2tt 2 T / c a b. The oscillat- 
ing part of the magnetization M = —dil/dB con- 
tains cosine functions of arguments 2irpF a /B = 
2irpf a x with p integer, and therefore the Fourier 
spectrum is made of the individual frequencies F a 
plus their harmonics. There is no mixing. In the 
CE, we will show that the oscillating part of the 
magnetization M = —dF/dB, where F = SI + 
n e fi is the free energy, contains the combinations 
fc/o + If i. The chemical potential at finite field 
is the implicit solution of the conservation equa- 
tion dil/dfi = —n e . Replacing the solution /x into 
Eq. (2) for ft, we can obtain the free energy, hence 
M. 

At zero temperature, we can solve the chemical 
potential equation exactly and obtain the exact to- 
tal energy E as a function of x. E can be expressed 
as E + b 2 G(x)/2, with G(x) an oscillatory func- 
tion, and E Q the zero field energy : 
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We used a combinatorial analysis to compute G. 
We had to find the [x] + 1 lowest LL between the 
2 bands, with the first [x] LL completely filled and 
the last partially. After tedious computations, we 
arrive at the following expression: 
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Eq. (4) is the main result for the two-band prob- 
lem at zero temperature. The oscillating part of 
the magnetization is simply m osc — n e G'(x)/2 = 
J2p A{F) sin(27rFx). In the double sum, it can 
happen that some terms diverge when I + l'a van- 
ishes in the denominator. But in this case, the 
quantity (I — l')a/(l + C\) = I is integer and the 
"interference" term 1 — cos(27rZ) in the numerator 
vanishes at the same time. These divergences are 
regularized by taking the finite limit of the quan- 
tity (1 — cos(e))/e 2 when e — > 0. 

Eq. (4) gives the amplitudes for the k th har- 
monics of the combinations F = l fo + hfi- h is 
a positive integer and Iq can be negative. These 
include "borbidden" frequencies like fi — fo- The 
amplitudes depend on whether fo / fi is equal to an 
irreducible ratio of 2 integers pi/po- In this case, 
we can find many combinations V fo + Ifi that are 
equal to kF. These occur when I = kl\ + mpi and 
/' = klo — mpo, with m positive integer. We then 
obtain 



A(kF) : 
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In the case where fo / fi is not rational, only the 
term m = has to be taken. We will only consider 
in the rest of the paper the case when fo/ fi is 
irrational, for simplication, and therefore only the 
first term m — on the right hand side of Eq. (5) 
contribute. Eq. (5) holds for any frequency except 
for the special case F = k(fo + fi)- 
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We obtain a simple result: the k th harmonics 
amplitude of fo + fi is the same as it would be 
in GCE case for a single band with fundamental 



frequency .fo + fi (except for a (— l) fc sign). The 
amplitudes for harmonics of f a are: 
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The difference with the GCE is that harmon- 
ics amplitudes of any frequency but fo + fi fall 
like 1/fc 3 instead 1/k. We conclude that the jumps 
in magnetization only come from the contribution 
of fo + fi, since its k th harmonics amplitude de- 
creases like 1/k. The amplitude of the "forbidden" 
frequency fi - /o is: 
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Figure 1 shows the numerical comparison 
between GCE and CE amplitudes. Analytical 
approximation based on thermodynamical re- 
lations [8] shows that the magnetization am- 
plitudes A(F) for combined frequencies F = 
hfo + hfi at finite temperature are proportional 
to Ro.i Q {T)R Ul {T)F/loh. Basically this is the 
product of GCE amplitudes for the individual fre- 
quencies lofo and hfi- For F = kf a , the authors 
found A(F) oc R a .kf a /k. At zero temperature 
the k th harmonics amplitudes of these frequen- 
cies decrease to 1/k in any case, which is not the 
exact case since we showed the dependence in k 
is 1/fc 3 . Moreover, the effect of commensurability 
between fundamental frequencies is not included. 
The finite low temperature dependence of the am- 
plitudes is therefore not understood yet. However 
their approximation may work beyond a crossover 
temperature below which corrective terms are 
needed. Indeed the chemical potential fluctuations 
decrease with increasing temperature and we ex- 
pect GCE and CE temperature reduction factors 
be equivalent for frequencies kf a . 

In conclusion, we computed the exact Fourier 
magnetization for the 2 band model at zero tem- 
perature. The expressions found are different from 
LK theory except for the harmonics of fo + fi- The 
k th harmonics amplitude of other frequencies de- 
creases like 1/fc 3 , and includes an interference term 
depending on the mass ratio. These results may be 
useful to check possible low temperature theories. 
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Fig. 1. Comparison between the absolute amplitudes 
in the CE (grey line) and GCE (black line) for 
ci = m*/m* = 1/2, A = 0.1, Ai = 0, n e = 1, /<, = 4/15 
and /i = 11/15. 
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